Discrete breathers in dc biased Josephson-junction arrays 
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We propose a method to excite and detect a rotor localized mode (rotobreather) in a Josephson- 
junction array biased by dc currents. In our numerical studies of the dynamics we have used 
experimentally realizable parameters and included self-inductances. We have uncovered two families 
of rotobreathers. Both types are stable under thermal fluctuations and exist for a broad range of 
array parameters and sizes including arrays as small as a single plaquette. We suggest a single 
Josephson-junction plaquette as an ideal system to experimentally investigate these solutions. 
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The phenomenon of intrinsic localization (intrinsic lo- 
calized modes or discrete breathers (DB)) is a recent dis- 
covery in the subject of nonlinear dynamics jjj. DB are 
solutions to the dynamics of discrete extended systems 
for which energy is exponentially localized in space. They 
appear either as oscillator localized modes, for which a 
localized group of oscillators librate; or rotor localized 
modes or rotobreathers, for which a group of oscillators 
rotate while the others librate ||. Recently, it has been 
found that DB are not restricted to periodic solutions but 
can also include more complex (chaotic) dynamics 

DB have been proven to be generic solutions in hamil- 
tonian |f and dissipative jf| nonlinear lattices. It is be- 
lieved that they might play an important role in the dy- 
namics of a large number of systems, such as coupled 
nonlinear oscillators or rotors. 

Though intrinsic localized modes have been the object 
of great theoretical and numerical attention in the last 
10 years, they have yet to be generated and detected 
in an experiment. Thus, finding the best system and 
method for the generation, detection and study of an 
intrinsic localized mode in a Condensed Matter system 
has become an important challenge Jg|.[7| . 

Josephson-junctions arrays are excellent experimental 
systems for studying nonlinear dynamics || . In this pa- 
per we propose an experiment to detect a rotating local- 
ized mode in JJ anisotropic ladder arrays biased by dc 
external currents |J. For this, we have done numerical 
simulations of the dynamics of an open ladder including 
induced fields [jl0| at experimentally accessible values of 
the parameters of the array. We also propose a method 
for exciting a rotobreather in the array. We distinguish 
between two families of solutions which present different 
voltage patterns in the array. Both types are robust to 
random fluctuations and exist over a range of parame- 
ter values and array sizes. Unexpectedly, we have found 
that many of the rotobreather solutions do not satisfy 
the up-down symmetry usually assumed for most types 
of dynamical solutions in the ladder. We also show that 
a DB solution can be most readily studied in a single 
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FIG. 1. Anisotropic ladder array with uniform current in- 
jection. Vertical junctions have critical current I cv and hori- 
zontal junctions I c h- 



plaquette. 

According to the RCSJ model, a Josephson junction 
is characterized by its critical current I Cl normal state 
resistance R n , and capacitance C. The junction voltage 
v is related to the gauge-invariant phase difference ip as 



<&o dip 
2ir~dt 1 



(1) 



where <£>o is the flux quantum. After standard rescaling 
of the time by r = \J &oC/2irI c , the normalized current 
through the junction is 



tp + Tip + sin (p, 



(2) 



here T represents a damping and is directly related 
to the Stewart-McCumber parameter (3 C = T~ 2 = 

2irI c CR 2 J<P - 

Our anisotropic JJ ladders (see Fig. ^) contain junc- 
tions of two different critical currents: I c h for the horizon- 
tal junctions and I cv for the vertical ones. Anisotropic 
arrays are easily fabricated by varying the area of the 
junctions. In the case of unshunted junctions, the crit- 
ical current and capacitance are proportional to this 
area. Due to the constant I c R n product, the normal 
state resistance is inversely proportional to the junction 
area. The anisotropy parameter h can then be defined as 
h = I c h/ lev = Ch/C v = R v /Rh- 

To write the governing equations of an anisotropic JJ 
ladder array with N cells, Fig. [l], we need to apply current 



1 



conservation at each node and flux quantization at each 
mesh. We are including self-induced magnetic fields so 
that flux quantization at mesh j yields 



(V Xip)j = -2nfj. 



(3) 



Here (V x <p)j = ip* + ipj + i — ip'j—ipj and it represents the 
circulation of gauge-invariant phase differences in mesh 
j = 1 through N. The self-induced flux through mesh j 
normalized by $o is given by fj. The resulting equation 
can be written compactly as, 



h{(p) + Tip* + simp*) = -A(V x ip) 

+ Tip] + sin$ = A[(V x ip) 3 - (V x ip)^} 
h{(p) + Tip) + sin <p)) = A(V x ip) h 



(4) 



where the open boundaries can be imposed by setting 
(V x ip) = (V x <p)n+i = in Eqs. ||. The system 
has four independent parameters: h, T, the penetra- 
tion depth A = &o/2irI cv L where L is the mesh self- 
inductance, and the normalized external current /. On 
writing Eqs.(^|) and (JJ) we assume zero external field 
and normalize by I cv . Non-zero external fields can be 
included in the model replacing the (V x (p)j terms in 
Eqs.(|) and (§) by (V x ip)j + 2 7 r/f :t where /? xt is the 
flux due to an applied external field, measured in terms 
of $ - 

The parameter values we will consider are based on 
Nb-Al2 x -Nb junctions with a critical current density of 
1000 A/cm 2 . Typical values of the Stewart-McCumber 
parameter and the penetration depth for arrays with h = 
1/4 are j3 c ~ 30 and A ~ 0.02. For the purposes of this 
work we will let T = 0.2, A = 0.02 and N = 8. 

Consider the h — limit in Eqs. |i| In this limit the 
vertical junctions behave as uncoupled damped pendula 
driven by an external current /. There, we can think of a 
configuration in which one or a few of the phases are ro- 
tating or oscillating around their equilibrium points while 
the others remain at rest. Thus rotor and/or oscillator lo- 
calized modes appear as solutions to the dynamics when 
the array is biased either by dc or ac external currents. 

For a single underdamped junction driven by a con- 
stant external current, the response measured in terms 
of dc voltage presents a hysteresis loop between the de- 
pinning and the retrapping currents. In this range the 
pinned (V — 0) and rotating (V ^ 0) solutions coex- 
ist. Then the rotobreather solution in the h — limit 
corresponds to a solution in which the phase of one of 
the vertical junctions is rotating while the other vertical 
junctions are at rest. 

As h is increased from zero the non-convex character 
of the coupling allows for the continued existence of ro- 
tobreathers in the system. Since a solution with a time 
increasing field cannot physically exist, the flux quantiza- 
tion condition (Eq. ||) implies that each cell with a rotat- 
ing junction must have at least one other junction which 
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FIG. 2. Schematic picture of rotobreather solutions: type 
A rotobreather in the ladder array (a), and plaquette (b); 
type B rotobreather in the ladder (c) and plaquette (d). Ar- 
rows are associated with rotating junctions and labels in (a) 
corresponds to graphs in Fig. || 

is rotating. Thus, for the single rotobreather solution 
one of the vertical and some of the horizontal neighbor- 
ing junctions rotate. Fig. |^ shows schematically simple 
rotating localized modes in a ladder and in the single pla- 
quette. These DB are amenable to simple experimental 
detection when measuring the average voltage through 
different junctions. 

Although the rotobreather solution can theoretically 
be continued from its h = limit by varying h, we have 
developed a simple method of exciting it in an array. 
This method should be experimentally reproducible and 
has three steps: (i) bias all the array up to the operating 
point (/ = I*); (ii) increase the current injected into 
one of the junctions to a value of the current above the 
junction critical current (I = I* + I > 1); (iii) go back 
to the operating point by decreasing to zero the value of 
this extra current /. Typical values of I* and / in our 
simulations are 0.6. 

We have checked the robustness of this method under 
fluctuations by simulating the equations of the ladder 
while adding a noise current to the junctions (this is the 
standard manner of including thermal effects in the sys- 
tem |ll]]). Thus we are able to excite DB in the ladder 
at some values of the parameters of the system. The so- 
lution showed in Fig. || was excited using this procedure. 

Henceforth, we are going to consider ladders with an 
even number of cells for which one vertical junction (the 
central one) is rotating. We will relabel this junction as 
3 = 0. 

Figs. [| shows a solution of a stable rotobreather in a 
JJ ladder. We plot the phase portrait (y*,^) of some 
of the superconducting gauge invariant phase differences 
of the array. The corresponding junctions are shown in 
Fig. ||(a). For clarity we have reduced the values of the 
phases to the (— tt, it] interval. We see that at this value 
of the penetration depth the solution is highly localized 
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FIG. 3. Example of a rotobreather in an 8 cell JJ ladder at 
h = 0.25, A = 0.02; T = 0.2 and I dc = 0.6. The only vertical 
rotating junction is the central one, j = 0, and we show the 
phase portrait of six of the phases of the array: (a) Lpg, (b) ip\, 
(c) ipo , (d) ip\, (e) ipo and (f) tp\ [Plotted phases are labeled 
in Fig. [](a)]. Also, in this case tpj = ipLj, = ~ V^-i- 

while three of the junctions describe a nearly sinusoidal 
rotation all the others oscillate with decreasing ampli- 
tudes. The average voltage through the three rotating 
junctions in the array is different from zero and equals 
to zero for all the other junctions. Fig. || shows the av- 
erage value of the induced field of the cells of the array. 
It decreases exponentially as fj ~ e _J ' ' 26 {j > and 

f-S = -fj-i) 

There are some surprising characteristics in this solu- 
tion. Current conservation in the open ladder implies 
ij = —ij. From Eqs. || and || we can see that <£>*■ = — (fj 
is a simple solution of the dynamics of the array and 
it corresponds to the up-down symmetry of the phases. 
All the previous theoretical approaches to the dynamics 
of the array (which include whirling modes, resonances, 
row-switching, etc.) and many of the numerical ones fo- 
cus on solutions which satisfy this up-down symmetry. 
However, looking at Fig. ||| we see that the rotobreather 
solution shown there does not comply with this simple 
symmetry; that is, </?*■ ^ —<pk although i* = — i b y 

We will distinguish between two families of single ro- 
tobreather solutions in the ladder which present different 
voltage patterns. The first family rotobreather A [see 
Fig.H(a)], is characterized by one vertical and two hori- 
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FIG. 4. The average value of the induced flux 

(f-j = —fj-i) at all the cells of the ladder for the roto- 
breather shown in Fig. ^| This average value decreases expo- 
nentially which is characteristic of DB solutions. 

zontal rotating junctions. Type A solutions have two pos- 
sible configurations. The two rotating horizontal junc- 
tions can be both in the same side, either top or bottom, 
as in Fig. ||(a); or one in the top and the other in the 
bottom. The second family, rotobreather B [FigJ|(c)], is 
characterized by one vertical and four horizontal rotat- 
ing junctions. The solution shown in Fig. || and ^ is a 
type A rotobreather. Up-down symmetric solutions be- 
long to family B but not all family B solutions satisfy 
this symmetry. 

Figs. |^ and || show a solution for which the scale of 
localization is smaller than one cell. Thus, it is natural 
to study the DB solution in the simplest ladder array, 
the single plaquette. Obviously, the concept of exponen- 
tial spatial localization is not applicable to the plaquette, 
but all the other characteristics of the solution remain. 
In particular we will also distinguish between type A and 
type B rotobreather solutions in the plaquette, which in 
this case correspond to one vertical and one horizontal 
rotating junctions [Fig.|^(b)], and one vertical and both 
horizontal rotating junctions [Fig.||(d)] respectively. The 
single plaquette biased by dc external currents, is then 
proposed as the simplest and most convenient experimen- 
tal system for detecting a rotating localized mode. The 
method for exciting the mode is also applicable to this 
system. 

An important experimental issue then becomes finding 
the region of existence of these DB solutions with respect 
to the system parameters in order to investigate the fea- 
sibility of designing an array in to detect a DB. To design 
an array we need to calculate the junction areas, so that 
the anisotropy needs to be known. Since different val- 
ues of the anisotropy affect the cell geometry, they also 
change the value of A. On the other hand, T is determined 
by the current density of the junctions and therefore it is 
fixed and independent of the geometry while the applied 
current can be easily changed while measuring. In order 
to make an optimal design we will fix the value of / and 
r to 0.6 and 0.2 respectively and study DB solutions in 



3 




FIG. 5. Maximum values of anisotropy for the existence of 
type A (solid circles) and type B (open circles) rotobreather 
solutions at different values of A in an 8 cell ladder (a) and a 
single plaquette array (b). F and I are equal to 0.2 and 0.6 
respectively. Lines serve as a guide to the eyes. 



the (h-X) plane of parameters. 

Type A and type B rotobreathers exist close to the 
h = limit. We then calculate the maximum value of 
the anisotropy for which a DB exists as a stable solu- 
tion of the dynamical equations for different values of A. 
Fig. H (a) and (b) show the result for the ladder and the 
single plaquette. The data were calculated by integrat- 
ing the equations of motion for the corresponding system 
with a small quantity of noise. We start with a type B 
rotobreather and h ~ 0.001. As we increase h, type B 
solutions become unstable and the solution evolves to a 
type A rotobreather. As we further increase h this roto- 
breather becomes unstable and the system usually jumps 
to either a pinned or a whirling state. To verify that our 
method is accurate, we have calculated Floquet multipli- 
ers for periodic rotobreather solutions and found results 
consistent with those shown in the Fig. |5[ 

We note that when doing this existence analysis of the 
solutions we find many different single-breather solutions. 
Most are periodic with different periods and amplitudes 
but there are some that appear to be chaotic [specially 
close to the A = 0.2 region in Fig. ||(a)]. A detailed study 
of the different bifurcations which include period dou- 
bling bifurcations to chaos is in current progress. There 
also exists a large family of different multibreather solu- 
tions, each one with its own domain of existence. 

Fig. H shows that at T — 0.2 and I = 0.6 type A solu- 



tions exist at larger values of the anisotropy than type B 
solutions. Also a simple inspection reveals a strong sim- 
ilarity between Figs. ||(a) and ||(b). The similarity can 
be easily understood. The rotobreather solution shown 
in Figs. ^ and |^ presents a mirror symmetry with re- 
spect to the rotating vertical junction: <pV = ip°_A, and 



<p_) = — vt-i' hi the case of solutions satisfying a mir- 



ror symmetry it is possible to map the dynamics of a JJ 
ladder for which the rotating junction is the central one, 
to the dynamics of a smaller JJ ladder for which the ro- 
tating junction is on one of the ends. Then, due to the 
localized nature of the DB solution the dynamics can be 
approximated by studying a single plaquette. When do- 
ing these transformations we need to rescale two of the 
parameters of the equations. Thus, results for the DB 
solution studied above present some similarities with the 
dynamics of a DB in a single plaquette when h p = 2hi 
and A p = 2A/. 

By establishing a criteria for the design of simple ex- 
periments to detect these intrinsic localized modes we 
hope to stimulate experimental investigations. 
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